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Abstract. A sample of nine bright Seyfert 1 and NELG type galaxies, observed with BeppoSAX, is analyzed to 
assess on a truly broad band basis (0.1-200 keV) the issue of the spectral contributions of Compton reflection and 
iron line fluorescence from circumnuclear gas. The empirical description adopted for the direct continuum is the 
commonly used power law with an exponential cut-off. The most direct test of the theoretical predictions, namely 
that the equivalent width of the line, W a , and the strength R of the reflection relative to the direct continuum 
are closely related to each other, gives a substantially positive result, that is their mean ratio is very close to 
expectation, and only a modest spread in the iron abundance seems implied. The existence of a steep correlation 
between R and the slope F of the power law is not confirmed. A weak evidence is found that the existence of a 
very shallow trend to increase on average with V cannot be altogether excluded in both R and W a , but needs to 
be tested with a larger sample. The energy Ef in the exponential cut-off spans a range from about 80 to more 
than 300 keV. A possible correlation is found, with Et increasing on average with V: if ignored, for instance by 
keeping Et at a fixed value in a sample study, it could be cause of artificial steepening in a correlation between R 
and T. 
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1. Introduction 

Observations with the Ginga satellite in the 2-18 keV 
band of a sample of bright Active Galactic Nuclei (Nandra 
& Pounds 1994), following the discoveries with the same 
satellite by Piro et al. (1990), Matsuoka et al. (1990) and 
Pounds et al. (1990), definitely confirmed the presence of 
reprocessed, in addition to the direct, power law emission, 
as predicted by e.g. Guilbert & Rees (1988), Lightman & 
White (1988), George & Fabian (1991), Matt et al. (1991). 
The reprocessed radiation consists of two main compo- 
nents. The first is seen as a broad excess rising above 
about 10 keV, due to Compton reflection of the power 
law by optically thick, low ionization matter, the second 
is K Q fluorescence of mainly "cold" iron at 6.4 keV. Both 
appear to be common features in Seyfert galaxies of type 
1 and intermediate between type 1 and 2 (objects of the 
latter type are often referred to as Narrow Emission Line 
Galaxies, NELG, their X-ray spectra are typically affected 
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by a substantial, order 10 22 23 H atoms cm 2 , neutral ab- 
sorbing column). 

Once the inclination angle i between line of sight and 
material is adopted, the amplitude of the reflection rel- 
ative to the power law can be converted into a measure 
R = fi/27r of the solid angle subtended by the reflector. 
Through a reanalysis of Ginga data for 23 radio-quiet 
Seyfert 1 and NELG, Zdziarski et al. (1999) found evi- 
dence of a strong correlation between R (estimated for a 
fiducial value of i — 30°) and the photon slope T of the 
power law, which they tentatively describe with a function 
where the average R (xT v , v about 12: R rises from about 
0.1 at r=1.7 to about 1 at r=2.1. They argue that the 
correlation is not an orientation effect, and, in the frame- 
work of the models where the production of the power law 
is associated with a very hot medium that comptonizes 
soft photons, they tentatively attribute it to a feedback 
of the additional photons from the reflection on the phys- 
ical state of that medium. At present, however, its real 
significance, at least in its quantitative terms, is in doubt, 
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because the existence of a strong correlation between these 
two variables in the fit of single observations is likely to 
bias substantially the result. Vaughan & Edelson (2001), 
for instance, on the basis of extensive simulations of data 
obtained with the PCA instrument onboard the RossiXTE 
satellite in an energy band similar to that of the Ginga 
data, have shown that a spurious correlation can be pro- 
duced. Done & Nayakshin (2001) proposed yet an alter- 
native possibility, namely that the correlation could be an 
artifact of the single-zone ionization (neutral in practice) 
model adopted for the reflector: however the simulations 
of Ginga spectra with their sophisticated model show that 
the correlation should be much shallower than the analyt- 
ical description quoted above. 

The equivalent width W a of the iron fluorescence line, 
which inevitably arises in the reflecting material, is ex- 
pected to be proportional to R. However, a line with W a 
up to about 100 eV, that is a substantial fraction of the ob- 
served values, in an optimal geometrical arrangement can 
be produced also by marginally thin gas (e.g. Wozniak et 
al. 1998). Thanks to the superior energy resolution, the 
observations with the ASCA satellite of a representative 
sample of objects (Nandra et al. 1997) have shown that 
in most cases the line is resolved and broader than the 
velocity spread of the optical Broad Line Region (BLR) . 
Moreover the first ASCA observation of MCG-6-30-15 
(Tanaka et al. 1995) revealed a profile consistent with the 
special and general relativistic effects predicted for a line 
originating in the innermost regions of a flat accretion disk 
(Fabian et al. 1989). Although the line width in this object 
turned out to be rather extreme compared to the others 
in the Nandra et al. (1997) sample, it is widely believed 
(e.g. Fabian et al. 2000) that the fluorescent gas is most 
likely associated with the optically thick accretion disk, 
and therefore the broad line W a should correlate with R. 
Unfortunately the ASCA band, being restricted to less 
than 10 keV, is inadequate to estimate R, and because of 
the variability which typically affects the continuum emis- 
sion of these objects (in addition to reported cases of vari- 
ations in the line, e.g. Weaver et al. 2001 and references in 
Fabian et al. 2000), the comparison of W a , as measured 
with ASCA, and of R as measured with Ginga, would lack 
significance. An alternative approach chosen by Lubinski 
& Zdziarski (2001) aims at looking for a correlation be- 
tween the broad line W a and T in ASCA observations of 
25 objects, which should correspond to the one between 
R and T. Their Fig. 1 shows indeed a trend, on average, 
in W a (single gaussian fit) to increase with T. This trend 
is quantified by fitting a relativistic diskline profile and an 
additional narrow gaussian profile to three average spec- 
tra, each one grouped according to the slope: the hardest 
group (r av =1.62) has broad W a —72^\l eV, the softest 
group (r av =1.94) has broad W a =124±26, with an average 
inclination angle (a product of the diskline fit) practically 
the same for the three groups, i ~ 45°. The additional 
gaussian has an average narrow W a of about 50 eV in 
all three groups. The authors interpret this component as 
originating in marginally thin gas: this is a possibility, as 



mentioned above, but the alternative one that it might 
instead be associated with highly reflecting material well 
beyond the accretion disk looks at present at least equally 
viable. For instance, the case of NGC 4051 observed with 
BcppoSAX in a state of exceptionally low intensity of its 
power law (Guainazzi et al. 1998) shows that the narrow 
line can be accompanied by the right amount of reflection, 
likely associated with a thick molecular torus well beyond 
the accretion disk, as predicted (e.g. Ghisellini et al. 1994) 
in Seyfert type 1 and 2 unified models (Antonucci 1993). 
The results in Lubinski & Zdziarski (2001) have been re- 
cently criticized by Yaqoob et al. (2002), and should be re- 
garded with some caution. Besides this criticism, we note 
that the trend in their Fig. 1 , mentioned above and which 
is based on individual measurements of W a and L, appears 
to be shallower than the one between R and T claimed by 
Zdziarski et al. (1999). 

In this context, the goal of this paper is to pro- 
vide an assessement of the observational contribution 
of the italian-dutch satellite BcppoSAX (Boclla et al. 
1997a). This satellite is characterized by: a very broad 
energy band, 0.1-200 keV, covered by its Narrow Field 
Instruments (NFI), namely the LECS (0.1-10 keV, 
Parmar et al. 1997), the MECS (1.8-10.5 keV, Boella 
et al. 1997b) and the PDS (13-200 keV, Frontera et al. 
1997); the energy resolution of the imaging ones (LECS 
and MECS) AE/E ~ 0.08(^/6 kcV)~i, which is interme- 
diate between that of ASCA and of Ginga at the iron line; 
the very effective control of the background in the PDS. 
These properties have shown their value for the estimate 
in several Seyfert 1 galaxies of the spectral parameters of 
interest for this paper, provided their flux in the 2-10 keV 
band is at least (2-3) xlO -11 erg cm -2 s -1 (e.g. Guainazzi 
et al 1999a, Perola et al. 1999). 

Thanks to the signal detected in the PDS band, the 
strength of the reflection is basically measured in the fit 
by the excess counts in the 20-50 keV, where the broad 
maximum of this component is located, over the power 
law whose slope is constrained by the spectral counts in 
the LECS and MECS. Determining the slope of the power 
law is therefore the main issue, and the confidence beyond 
the purely statistical errors is related to the confidence 
that can be associated, object by object, to the treatment 
of additional spectral components, such as the cold and 
warm (ionized gas) absorbers and a generic soft excess. 
In addition, the power law itself may be affected by a 
high energy "break" (as first seen in NGC 4151, Jourdain 
et al. 1992), which we shall, as it is often done in the 
literature, empirically describe with an exponential factor, 
exp(— E/Ef ). Since the relative strength of the reflection 
below and at its broad maximum is very sensitive to the 
Compton recoil in the scattering of higher energy photons, 
hence to E{ , the ability to estimate both R and Ef thanks 
to the wide band and the efficient background control of 
the PDS (collimator used in rocking mode) is an asset 
of BcppoSAX. Nonetheless caution must be retained on 
single observation measurements of the three parameters 
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r, R and Ef, because systematic errors on V, arising from 
an imperfect description of the spectrum at low energies, 
could be substantial compared to the purely statistical 
errors. 

When the iron line is resolved, it might seem more ap- 
propriate to use models, both for the line and the associ- 
ated reflection, which incorporates Doppler and relativis- 
tic effects. This approach would be physically more self- 
consistent, and in principle provide in addition the incli- 
nation i. However, as confirmed by observations with the 
Chandra and XMM-Newton satellites (detailed in Sect. 3 
for our objects), the presence of a substantial narrow con- 
tribution is a rather common property, which is difficult to 
separate from the broad contribution with the BcppoSAX 
resolution, thus making this approach very tricky, particu- 
larly for the objects in our sample, where the line width is 
modest, a < 0.4 keV. We have therefore chosen to follow 
a conservative approach, where the model adopted for the 
reflection component is static, and the line is described 
as a single gaussian. We cared, though, to verify the out- 
come of adding to the model a narrow line at 6.4 keV: in 
general no improvement in the fit is achieved, moreover, 
when both a narrow and a broad component are found, 
their combined W a does not differ significantly from that 
obtained with a single gaussian. 

Despite their BeppoSAX observations fulfilled the se- 
lection criteria given in Sect. 2, in order to limit the inci- 
dence in our sample of the systematic errors mentioned 
above, we prudently excluded objects with a too com- 
plex and variable absorber, like NGC 4151 (Piro et al. 
2001), NGC 3516 (Costantini et al. 2000, Guainazzi et al. 
2001) and NGC 2992 (Gilli et al. 2000). The observation 
of MCG-6-30-15 (Guainazzi et al. 1999b) confirmed the 
ASCA finding (Tanaka et al. 1995) of a very broad and 
asimmetrical iron line, and was also excluded because it 
could not be treated homogeneously with the rest of the 
sample. 

The paper is organized in the following manner. Sect. 
2 describes the sample of the observations used. Sect. 3 
describes the spectral fitting procedures and the results. 
Sect. 4 illustrates the distribution of the parameters in 
connection with the reprocessing. Sect. 5 is devoted to 
discussion and conclusions. 



2. The sample data 

The sample contains observations of objects with flux in 
the 2-10 keV band greater than 3xl0 -11 erg cm~ 2 s _1 , 
duration longer than 50000 s in the MECS. For some ob- 
jects, with separate observations shorter than 50000 s, af- 
ter verifying with a hardness ratio analysis, as well as in- 
dependent fits, that no spectral changes had occurred of 
such entity as to impair the reliability of the outcome rele- 
vant to the present paper, we merged two (NGC 5506, Mrk 
509) or three observations (IC4329A(2)). Several of these 
observations are the subject of papers already published, 
whose reference will be given in Sect. 3. 



Table 1 lists for each object: the observation start- 
dates and the good exposure times in the LECS, MECS 
and PDS, the rcdshift and the galactic hydrogen column, 
N g , adopted, the flux in the 2-10 keV band as observed, 
F Q , corrected for any absorption as from the model fit, F na , 
and its corresponding luminosity (H Q ~75 km s _1 Mpc -1 ). 

The reduction procedure and screening criteria 
adopted are standard (see e.g. Guainazzi et al. 1999a). 
Concerning the PDS data, the Variable Rise Time (VRT) 
option was generally adopted, except for the observation 
of NGC 4593, where the Fixed Rise Time (FRT) was used. 
The PDS net spectra were obtained from direct subtrac- 
tion of the off- from the on-source products. We verified 
that for none of the sample objects there was a signifi- 
cant contamination by hard sources in both the on- and 
the off-source field of view. Concerning the data from the 
LECS and the MECS images, the background was sub- 
tracted using spectra from blank sky files at the source 
position in the detectors. 

The exposure times in the LECS are much shorter than 
in the MECS, because the former instrument is operated 
only during the night-time fraction of each orbit. Since 
the incident flux is generally recorded as variable during 
each observation, to take care of this inhomogeneity in the 
time coverage, in the fits the normalization of the LECS 
to the MECS was left as a free parameter. The normal- 
ization of the PDS to the MECS was instead held fixed 
at 0.8 (VRT selection) or 0.86 (FRT selection, Fiore et al. 
1999). For the LECS data we used the version Sept. 1997 
of the transfer matrix, and limited the range to 0.1-4.0 
keV, where this matrix is more reliable. The energy bins 
chosen correspond to about one third of the instrument 
resolution, which itself is energy dependent. 

3. Spectral analysis 

The analysis of the spectral counts was performed us- 
ing the software package XSPEC (Arnaud 1996, ver- 
sion 11.0.1). The Baseline Model Spectrum (BMS) in- 
cludes: a photon power law (PL) with exponential cut-off, 
AE~ V exp(—E/Ef), together with a reflection component 
(RC) from a cold slab isotropically illuminated by the PL 
photons, subtending the solid angle R = Q./2-n and with 
inclination cos i to the line of sight (Magdziarz & Zdziarski 
1995, module pexrav); a uniform cold absorber column 
at the source, N s , in addition to N g given in Table 1; a 
uniform warm absorber at the source, as described in the 
module ABSORI, with column N w and ionization parame- 
ter X\ — Li/nD 2 erg cm s _1 , where n is the gas density 
and D its distance from the source of ionizing luminosity 
Li in the interval 5 eV to 20 keV (for L;, assumed with a 
power law spectrum, we set the slope equal to T; the tem- 
perature of the absorbing gas was fixed as a rule at 3xl0 4 
K); a gaussian line to represent the iron K Q fluorescence, 
with energy E a , width a a , intensity I a , accompanied by a 
gaussian of identical width, to represent the iron Kp (with 
Ef} and Ip set equal to l.lE a and 0.11/q,, George & Fabian 
1991). For both reflection and absorption we adopted the 



4 G. C. Perola et al.: Compton reflection and iron fluorescence in Seyfert 1 galaxies 

Table 1. The sample 



Source 
name 


otart 
date 


t pre 

T oxp (s) 


11 If T?pC 

T oxp (s) 


Toxp (s) 


z 


N g 

(xlO 20 cm-2) 


F 

(2-10 keV) 


(2-10 keV) 


(2-10 keV) 


MCG 8-11-11 


1998-Dec-07 


29840 


74865 


72472 


0.020 


20. 3 1 


5.61 


5.73 


43.62 


MCG-5-23-16 


1998-Apr-24 


35825 


76941 


65048 


0.008 


8.4 2 


8.73 


10.00 


43.06 


NGC 3783 


1998-Jun-06 


69135 


153870 


140814 


0.009 


9.6 2 


6.44 


6.86 


43.00 


NGC 4593 


1997-Dcc-31 


31426 


84057 


76934 


0.009 


2.0 1 


3.72 


3.79 


42.74 


IC 4329A (1) 


1998-Jan-02 


25097 


81826 


75214 


0.016 


4.5 1 


13.15 


13.83 


43.81 


IC 4329A (2) 


1998-Jul-17 
1998-Jul-21 
1998-Jul-27 


38245 


106140 


94772 






10.68 


11.18 


43.71 


NGC 5506 


1997- Jan-30 

1998- Jan-14 


25101 


77491 


70842 


0.006 


4.3 2 


7.20 


9.57 


42.79 


NGC 5548 


1997-Aug-14 


85444 


312810 


239320 


0.017 


1.6 2 


3.55 


3.61 


43.27 


Mrk 509 


1998-May-18 
1998-Oct-ll 


41786 


87834 


82102 


0.035 


4.4 2 


5.66 


5.69 


44.10 


NGC 7469 


1999-Nov-23 


99114 


249780 


243020 


0.017 


4.8 1 


3.73 


3.76 


43.29 



a Redshift values from Veron-Cetty & Veron (1993) 

b Galactic H column from: (1) Elvis et al. (1989) and (2) Murphy et al. (1996) 
c F : Flux as observed; F na : corrected for absorption; both in units of 10~ n 
d Log of luminosity corrected for absorption in erg s _1 (H =75 km s _1 



erg cm 
Mpc" 1 ) 



element abundances in Anders & Grevesse (1993). All fits 
were performed with cos i fixed equal to 0.9 (i = 26°). 
Alternatively cos i was left free with R fixed equal to 1. 
Because the dependence on cos i of the RC shape is mod- 
est, all other parameters turned out identical in practice 
to those obtained with the angle fixed, and the differences 
in x 2 are insignificant. 

The results on the parameters of main interest for this 
paper are collected in Table 2. The line parameters are 
given in the source frame, and the line intensity I a is also 
expressed as an equivalent width, W a . The errors repre- 
sent the 90% confidence interval for two parameters of in- 
terest (Ax 2 =4.61, Lampton et al. 1976), a realistic choice 
dictated by the statistical correlation which exists between 
couples such as (r, R), (T, Ef), (a a , I a ). 

To be noted that in general the value of E a is perfectly 
consistent with 6.4 keV, the value for "cold" iron (in agree- 
ment with previous findings, e.g. Nandra & Pounds 1994, 
Nandra et al. 1997, Lubihski & Zdziarski 2001). In two ob- 
jects (NGC 5506 and Mrk 509) the (per se marginal) evi- 
dence of a larger value finds explanation in recent XMM- 
Newton results, reported below. 

Next we briefly describe the case of each source, in- 
cluding the best fit values of the other BMS parameters 
and additional spectral components adopted to improve 
the quality of the fit when the BMS gave an unaccept- 
able x 2 . Generally Xi is very poorly constrained, hence 



we give only its best fit value. Reference is also given to 
some recent results on the iron line. 

MCG 8-11-11. This observation was published by 
Perola et al. (2000), to which we refer for the choice of 
neglecting in the fit the LECS counts below 0.4 keV. 
N s is negligible compared to N gl the warm absorber has 
AT w =(6^)xl0 20 cm" 2 with -7 erg cm s _1 . The iron 
line is not resolved, the upper limit on its width is com- 
parable to that of the lines resolved in our sample. 

MCG-5-23-16. The cold absorber in this NELG is sub- 
stantial, A s =(1.6±0.9)xl0 22 cm~ 2 ; no signature is de- 
tected of a warm absorber. In ASCA (Weaver et al. 1997) 
and RossiXTE (Weaver et al. 1998) observations (com- 
parable in flux to this one) a strong (W a ^250 eV) and 
broad iron line was consistently found. In the BeppoSAX 
observation the line is much narrower than a a ~ 0.4 keV, 
the value obtained from a single gaussian description in 
both the ASCA and the RossiXTE data sets (Weaver et 
al. 1998); the line intensity is about twice that of the nar- 
row gaussian, but about one third of the sum of the narrow 
and broad gaussians, as measured by Weaver et al. (1997) 
in their fit to the ASCA data. While this is probably yet 
another case of variations in the iron line (see also Smith 
& Done 1996 for a formally significant difference in W a 
between two Ginga observations made a few days apart), 
we note that, contrary to expectation, the amount of re- 
flection measured with RossiXTE (for a choice of £^=200 
keV, the closest to our best fit value of the three adopted; 
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Source 
name 


r 


E { 
(keV) 


R 

(cos i=0.9) 


E a 
(keV) 


(keV) 


1° 


(eV) 


X 2 /dof 


cos i 
(R=l) 


X 2 /dof 


MCG 8-11-11 


, oc+0.09 
i - SD -0.05 


166l?. 4 5 


1 09+ 80 


6.49±0.14 


< 0.37 


6.9i 3 2 -; 4 


117+ 63 


132.3/141 


o.96ir 4 4 


132.6/141 


MCG-5-23-16 


1.81±0.05 


1471™ 




6.43±0.10 


< 0.28 


lO.ol 3 , 9 


96l 3 2 7 


139.5/135 


0.511"; 24 


140.7/135 


NGC 3783 


1.77±0.04 


156lf 


0.631°;^ 


6.38±0.09 


97+ 15 
U - £l -0.13 


11 7+2-8 
"■'-2.7 


16111? 


150.4/148 


46+° 19 


148.0/148 


NGC 4593 


1 04+ 006 

*■•"*— 0.05 


> 222 


i.26l ; 40 


6.41±0.16 


34+ 036 




23ll^ 2 1 


148.2/146 


> 0.63 


149.0/146 


IC 4329A (1) 


1.89±0.04 


dZO -105 


0-63t°;iI 


6.43±0.17 


„ on + 0.30 

u.au_ 26 




104+ 41 


146.9/141 


0.49±°;£ 


148.2/141 


IC 4329A (2) 


1.90±0.05 


262l^ 4 


7-?+ - 22 


6.54±0.14 






123±« 


154.4/140 


0.58l«;f 4 


153.9/140 


NGC 5506 


9 n9 + 0.09 


> 298 


1 20+ 45 




21+ 015 


14.21 4 ; 2 


152±« 


85.3/82 


> 0.75 


85.9/82 


NGC 5548 


1 62+ 04 
LO - i -0.05 


147l 6 3 4 3 


o-64±S:?S 


6.41±0.07 


< 0.13 


4-31- 


106lf 3 


175.1/148 


-?7+ - 21 
u - d ' -0.12 


176.7/148 


Mrk 509 


, r„ + 0.09 

1 - oo -0.08 


ul -20 


o-68±S:»S 


6 64 +0 34 


< 0.87 


5.5I2I 


no + 76 


135.7/132 


45 + - 50 

u ' d -0.33 


137.3/132 


NGC 7469 


1 OO + 0.05 
1 - 88 -0.07 


164+" 6 

— DO 




6 44+° 12 
-0.11 


< 0.47 


5.211- 


1391^ 


168.4/148 


U - db -0.20 


168.0/148 



a Line intensity in units of 10 5 photons cm 



- 2 s- 1 



T is identical to our best fit value; Weaver et al. 1998) 
is definitely not larger than in the BeppoSAX observa- 
tion. Finally we note that a Chandra HETGS observa- 
tion, briefly described in Weaver (2001), shows a narrow 
(FWHM less than 3000 Km s" 1 ) line with W a ~ 90 eV, 
similar to our value. 

NGC 3183. The BMS fit is poor (x 2 =208/150), with 
an excess in the residuals at low energies. With the ad- 
dition of a line the x 2 reduces to 160/147, with £q=0.59 
keV and Wi=140 eV: its strength is however almost ten 
times larger than that of the OVII emission lines detected 
by Kaspi et al. (2001) in a Chandra HETGS observation 
(when the source was fainter by only about 15%). A for- 
mally better improvement is achieved if, instead of a line, 
a black body is added, provided the temperature of the 
warm absorber is set equal to 10 5 K: the parameters given 
in Table 2 refer to this fit. N B is negligible compared to N s , 
iV w =(9.6l5:4) xl ° 21 cm" 2 , with X { -15 erg cm s _1 , the 
black body kT=0.l8t° ° o l keV. The iron line is definitely 
resolved. The Chandra HETGS observation (Kaspi et al. 
2001) has revealed a narrow line at 6.4 keV, with an upper 
limit on the width (a less than 0.030 keV) that places its 
origin beyond the BLR: its intensity of (6.6±2.1) x 10~ 5 
cm~ 2 s _1 is about 55% of that measured by us. If we in- 
clude in the fit a narrow, in addition to the broad line, no 
improvement is attained in \ 2 an d the intensity of the nar- 
row component is not well constrained, such that a value 
identical to that from the Chandra observation is accept- 
able. In a recent publication (De Rosa et al. 2002) on this 
BeppoSAX observation, the same model adopted by us 
(their Model E) yields somewhat different values, in par- 
ticular r=1.86±0.03, i?=0.71±^28, TT Q =210±45, together 
with N w about twice larger. The discrepancy is due to the 
different code adopted to describe the warm absorber, and 



it is not such as to influence the conclusions of the analysis 
following in Sect. 4. We note on the other hand the good 
agreement in our values of T and iV w with those obtained 
by Kaspi et al. (2001) in fitting the "line free zones" of 
the HETGS spectral data with a power law and a warm 
absorber. 

NGC 4593. This observation was published by 
Guainazzi et al. (1999a). N s is negligible compared to N g , 



and7V w =(2.3 



+1.1 



0.9) x 10 



21 



cm with Xi —9 erg cm s . The 



iron line is resolved, W a together with R are the largest 
in our sample. 

IC 4329A. Observation (1) was published by Perola 
et al. (1999). The average flux of observation (2) is 20% 
lower, yet there are no statistically significant differences 
between the two sets of parameters (the same holds for 
the three pointings merged, when fitted individually, de- 
spite a 25% difference between the lowest and the highest 
of their flux levels). The mean values of the absorbers 



^ 21 cm" 2 



are: JV B =(2.7±1.0) x 10 21 cm" 2 , AT W =(2.8I^) x 10 
with X[ ~5 erg cm s _1 . The iron line is resolved. Given the 
strength of the source (the brightest of the sample) and 
the correspondingly good statistics, Perola et al. (1999) 
attempted a fit with a narrow gaussian plus a relativis- 
tic disk profile, and obtained marginal evidence for the 
existence of the narrow component. Another marginal de- 
tection of the narrow component in a simultaneous ASCA 
and RossiXTE observation is presented by Done et al. 
(2000). In a XMM-Newton observation with the EPIC in- 
strument (Gondoin et al. 2001), at a flux level 25% higher 
than in observation (1), the narrow (a less than 0.06 keV) 
line is detected, with W Q =43±1 eV; a broad component 
is seemingly absent, but unfortunately the authors do 
not provide an upper limit on its equivalent width. A si- 
multaneous and relatively short (34000 s in the MECS) 
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BcppoSAX observation is used, together with the EPIC 
data, by Gondoin ct al. (2001) to estimate r=1.93±0.03 
and i?=l.l±0.3: we believe that the estimate of R is incor- 
rect, because obtained with a PDS to MECS normaliza- 
tion factor equal to an improbable value of 0.7; by fitting 
the BcppoSAX data, as retrieved from the public archive, 
with this factor equal to 0.86 (for FRT selection, see Sect. 
2), we obtain best fit values of i?=0.54 and r=1.91, fully 
consistent with those in Table 2. 

NGC 5506. The two observations here merged have an 
almost identical average flux level, and when fitted indi- 
vidually yield similar values of the spectral parameters. 
The BMS fit with all three instruments gives an unac- 
ceptable x 2 =180/117, due to an excess in the residuals 
well below the strong photoelectric cut-off. Since the pa- 
rameters of interest do not crucially depend on how this 
excess is modeled, we have chosen to present in Table 2 the 
results obtained excluding the LECS. This NELG has a 
substantial iV s =(3.7±0.2)xl0 22 cm" 2 , and no signatures 
of a warm absorber. The iron line is resolved and E a is 
marginally inconsistent with 6.4 keV. This line is likely 
a blend of the two components recently revealed in a si- 
multaneous XMM- Newton (with EPIC) and BeppoSAX 
observation (flux level similar to ours, Matt et al. 2001), 
one unresolved at 6.4 keV, the other resolved (a about 
0.25 keV) at 6.75 keV, the two with comparable intensi- 
ties whose sum is consistent with the one in our obser- 
vation. Matt et al. (2001) present arguments, centred on 
the interpretation of the iron edge seen at 7.1 keV, that a 
substantial fraction of the RC in this object is very likely 
associated with the narrow line component. Evidence sup- 
porting this interpretation can be found in a spectral vari- 
ability investigation conducted with RossiXTE by Lamer 
et al. (2000). 

NGC 5548. This observation was published by Nicastro 
et al. (2000). They find internally evidence of a corre- 
lation between T and the flux (previously discovered in 
Ginga observations by Magdziarz et al. 1998, confirmed 
with RossiXTE observations by Chiang et al. 2000), which 
has some impact on the quality of our fit. In a Chandra 
LETGS observation Kaastra et al. (2002) find evidence of 
a very complex warm absorber, of a soft excess and of low 
energy emission lines. Following Kaastra & Barr (1989), 
who first reported a soft excess in EXOSAT observations, 
they model it as a "modified black body" with tempera- 
ture kT —0.1 keV. They suggest that the NVII emission 
line at 0.5 keV could be identified with the low energy 
line detected by Nicastro et al. (2000) in the BeppoSAX 
observation. We therefore added to the BMS a soft com- 
ponent, as described by the module DISKBB in XSPEC, 
with kT fixed at 0.1 keV, and a line with E\ fixed at 
0.5 keV. The fit yields Wi=76±^ eV, compatible with 
the Chandra result. N s turns out negligible compared to 
N g , while 7V w =(2.4±0.7)xl0 21 cm" 2 (with X\ -14 erg 
cm s _1 ) is marginally consistent with the column den- 
sity of the dominant component in the warm absorber as 
estimated by Kaastra et al. (2002). The iron line is unre- 
solved and significantly narrower (less than 0.13 keV) than 



in any of the other sample objects. A Chandra HETGS 
observation (Yaqoob et al. 2001), with the source at a 
flux level about 30% lower, has revealed a narrow line 
at 6.4 keV with a a =41^l'l eV, corresponding to a veloc- 
ity spread consistent with an origin in the outer BLR (a 
substantial contribution to the line from thick matter be- 
yond that region is not excluded). The intensity of this 
line matches the one measured by us. Using the same pre- 
scriptions for a relativistic disk line used by Yaqoob et al. 
(2001) to obtain an upper limit of 7xl0~ 5 cm~ 2 s _1 , we 
find a more stringent limit of 4xl0~ 5 cm~ 2 s _1 : this is a 
factor about 2 less than in earlier ASCA records of a sub- 
stantially broader line, as reported in Yaqoob et al. (2001) 
or in Chiang et al. (2000). In a spectral variabilty study 
with RossiXTE, Chiang et al. (2000) found the intrigu- 
ingly paradoxical (for the reprocessing paradigm) result 
that W a is anticorrelated to R. 

Mrk 509. The two observations and the results from 
their merging were published by Perola et al. (2000): in 
that paper it is shown that the BMS fit leaves the line 
width practically unconstrained (about 3 keV) and that 
in this respect the situation changes radically if a soft com- 
ponent is added in the form of a power law (preferred to 
a black body). This component has r s =2.5^0 2 an d nor- 
malization 1.1 xl0~ 2 cm~ 2 s _1 keV -1 at 1 keV. N s is neg- 
ligible compared to N g , and jV w =(4jJ) x 10 20 cm~ 2 with 
X{ —16 erg cm s _1 . The line width is poorly constrained 
and the line energy is marginally inconsistent with 6.4 
keV. A XMM- Newton observation (with EPIC, Pounds et 
al. 2001) has revealed the presence of a resolved cold and 
of a broader ionized (6.9 keV) components. The W of the 
former is about twice that of the latter, their sum is consis- 
tent with the value measured by us. The source, however, 
was almost 3 times fainter than in our observation, thus 
the blend must have experienced a change in intensity by 
a comparable factor. 

NGC 7469. The BMS fit gives a marginally acceptable 
X 2 = 191/151, but, like in the case of Mrk 509 as illustrated 
in Perola et al. (2000), the line width is practically uncon- 
strained (about 2 keV). The addition of a soft excess in 
the form of a power law changes the situation radically. 
The slope of the soft power law is r s =3.3±0.5, its normal- 
ization 5.5xl0~ 3 cm~ 2 s _1 keV -1 at 1 keV. N s is negligi- 
ble compared to N g , and N w =(6.7±f%) xlO 20 cm" 2 with 
Xi —0.13 erg cm s _1 . The line is not resolved, the upper 
limit on its width is comparable to that of the lines re- 
solved in our sample. In a spectral variabilty study with 
RossiXTE, Nandra et al. (2000) detected changes in the 
line intensity which are not correlated with the absolute 
normalization of the reflection component, but they them- 
selves doubt the reliability of their estimate of the second 
of the two parameters. We note that George et al. (1998) 
did not find evidence of a soft excess in their analysis 
(above 0.6 keV) of a ASCA observation with a flux 25% 
lower than in our data: thus the additional soft compo- 
nent, that we suggest here to provide a better description 
of the continuum of this observation, should be variable. 
We further note that the value of T in Table 2 is in good 
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agreement with r=1.92±0.02 obtained by Nandra et al. 
(2000) from the integrated fit of the data in the study men- 
tioned above, with an average flux only 10% lower than in 
our observation. 

4. An analysis of the parameters 

This Section is devoted to illustrate the distribution of the 
parameters in the sample, in relation with the reprocessing 
of the PL. 



4.1. W a and R 

Irrespective of model details, under the hypothesis that 
the line emission is entirely associated with optically thick 
material, to a good approximation W a should correlate 
linearly with the value of R obtained for an arbitrarily 
fixed inclination angle. Apart from the one introduced by 
the measurements errors, an intrinsic scatter in W a about 
the correlation can be caused by a variance in the iron 
abundance (for differences within a factor 2 of the cosmic 
value, the dependence of W a on the iron abundance is 
approximately linear, see Matt et al. 1997). Furthermore, 
for a given R, W a may differ according to the value of 
T, because it is proportional to the ratio of the photon 
number above 7.1 keV responsible for the fluorescence to 
the photon flux density at the line energy. The plot in Fig. 
1 appears scatter dominated, and because of the rather 
large errors on both parameters little can be said about 
individual deviations from the hypothetical linear trend. 
The Spearman rank-order (Press et al. 1992) test, which 
does not take the errors into account, yields a correlation 
coefficient which is positive, r s =0.43, but with a large, 
0.2, significance level for this quantity deviating from zero. 
Following another approach, where the error bars are also 
ignored, but which seems more useful in this context, we 
can compute for each object the ratio of the best fit value 
of W a to the value expected, given the T of the PL and 
the best fit value of R for cos z=0.9. We did that by using 
the T and angular dependence of W a given in George & 
Fabian (1991), rescaled to the Anders & Grevesse (1993) 
abundances (see Matt et al. 1997). We obtain 0.97 for 
the mean value, and 0.29 for the standard deviation of 
the sample. A similar outcome (mean value 0.88, standard 
deviation 0.27) obtains when we proceed instead with R 
fixed equal to 1 and the best fit values of cos i to estimate 
the expected W a . In our sample, then, on average the 
observed W a is of the right amount expected from the 
simultaneously measured strength of the RC relative to 
the PL for the cosmic value of the element abundances. 
The standard deviation, being inevitably affected by the 
large errors on both W a and R can be considered a limit 
on the variance of the iron abundance in our sample. 

Concerning R, we note the concentration (6 out of 9 
objects) of best fit values in the interval 0.50-0.73, or the 
corresponding concentration of cos i (for R=l) in the in- 
terval 0.35-0.58. The absence of objects with cos i less 
than 0.35 (inclined more than about 70°) can be under- 
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Fig. 1. W a plotted versus R. 
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Fig. 2. R plotted versus T. 

stood in the frame of the Seyfert type 1 and type 2 unified 
models (Antonucci 1993), where the first type takes the 
appearence of the second at large inclinations. The paucity 
of objects with small, relative to intermediate inclinations 
could well be a statistical fluctuation, given the modest 
sample size, but it may reflect intrinsic differences in R, 
with objects where the intrinsic R is indeed substantially 
less than unity. 
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Fig. 3. W a plotted versus T. The line, arbitrarily normal- 
ized, describes the trend expected (derived from George 
& Fabian 1991). 



4.2. R and W a versus T 

In a diagram of R versus T, scatter is expected due 
to differences in orientation. Fig. 2 shows scatter in R 
only above r=1.8, where the largest values are confined 
together with values comparable to those found below. 
While it could be due to a statistical fluctuation in the 
sampling of the angles, this mere fact provides some evi- 
dence of a trend in R to increase on average with increas- 
ing r. This is reflected in the Spearman rank-order cor- 
relation coefficient r s , which is positive and equal to 0.67, 
with a significance level for r s deviating from zero equal to 
0.03. The marginal significance of this evidence is further 
weakened by the consideration (see Sect. 1) that the two 
variables R and T remain far from independent also in the 
fits to our very broad band data, as can be appreciated 
with a simple exercise: by fitting the data of each obser- 
vation with r fixed at the two extremes of its error bar, 
the values of R obtained differ typically by as much as 60 
to 80% of their 90% confidence intervals. In any case, the 
trend in our sample (admittedly much smaller, 9 objects 
as opposed to 23) is definitely much shallower than the 
correlation found by Zdziarski et al. (1999; see also Sect. 
4.3 for a supplementary cause likely to contribute to the 
difference) . 

In addition to the one due to differences in the ori- 
entation, a scatter due to differences in abundance is ex- 
pected in a graph of W a versus V. A weak anticorrelation 
is predicted, with W a decreasing with increasing T, as 
mentioned above. The line (arbitrary normalization) de- 
scribing this effect (deduced from George & Fabian 1991) 



is plotted in the graph of W a versus V in Fig. 3. Notably 
W a shows instead a weak tendency to increase, on aver- 
age, with r, which resembles that obtained by Lubinski 
& Zdziarski (2001, their Fig. 1). This is reflected in the 
sign of r s =0.57, which, contrary to the prediction, is pos- 
itive. The significance level for r s deviating from zero is 
however 0.08, three times larger than in the previous case. 
We note, though, that in the fits the variable W a is less 
dependent on T than the variable R: the same exercise 
performed for R, and just mentioned, yields differences 
in W a which are at most, and typically less than 35% of 
their 90% confidence interval. Thus we are inclined to re- 
gard the evidence of a shallow trend in W a increasing on 
average with T, albeit very marginal, no more so than it 
appears to be for R. 



4.3. Ef versus T 

The plot of Ef versus T in Fig. 4 appears to confirm the 
existence of a correlation, with Ef increasing on average 
with T, noted first in BcppoSAX data by Piro (1999) and 
then by Petrucci et al. (2001, their Fig. 6). Although the 
variable Ef in the fits, very much like R, is rather strongly 
dependent on T, and therefore one should not disregard 
the possibility that our result might be biased by this ef- 
fect, we note that in this case the correlation coefficient 
r s is equal to 0.88, and the significance level for r s deviat- 
ing from zero is comfortably lower and equal to 0.0007. 
Discussing this correlation in terms of the mechanisms 
that give rise to the main spectral component goes beyond 
the scope of this paper (but see Petrucci et al. 2001). Our 
purpose here is to point out that, so long as one adopts 
the empirical representation of the direct continuum as a 
power law with an exponential cut-off, the existence of 
this type of correlation could, if ignored, artificially con- 
tribute to a correlation between R and T: in fact, as Ginga 
data do not allow to evaluate Ef, by either neglecting the 
exponential cut-off or keeping it at a fixed value irrespec- 
tive of T (Zdziarski et al. 1999 choose Ef = 400 keV fixed), 
one would systematically obtain lower values of R for lower 
values of T, for the reason already mentioned in Sect. 1. 
As an example of the interplay between Ef and R, we 
show the case of NGC 3783, whose T has an intermediate 
value in our sample: after excluding the cut-off we obtain 
X 2 =212. 1/149 (A% 2 =61.7 demonstrates the high signifi- 
cance of the cut-off in the fit given in Table 2) and a best 
fit value of i?=0.38 as opposed to _R=0.63. Although the 
effect is sizeable, it alone is probably insufficient to ex- 
plain the difference in steepness of the correlation found 
in the Ginga sample by Zdziarski et al. (1999) and in our 
sample, but it does represent an explanation of the offset 
between Ginga and BeppoSAX results noted, and tenta- 
tively attributed to calibration differences, by the same 
authors. 
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Fig. 4. The cut-off energy Ef plotted versus T. 
5. Discussion and conclusions 

In confronting theory of reprocessing in AGN and obser- 
vations, first comes a test of the predicted correlation be- 
tween the strength of the reflection and that of the iron 
line fluorescence relative to the primary photons in the 
PL, from material mainly al low ionization, as indicated 
by the observed energy of the line. Our test, while improv- 
ing on it, does confirm the result originally obtained by 
Nandra & Pounds (1994), who compared the mean sam- 
ple values of R and W a with expectations from a face-on 
slab. We take properly into account the dependence of W a 
on T in each observation and show that, be R intrinsically 
equal to one or otherwise, the amount of fluorescence is 
close to what expected from the strength of the reflec- 
tion, and that a modest variance in the iron abundance 
is likely present. In doing that, we assume that the fairly 
substantial "narrow" contribution to the line, which from 
Chandra and XMM-Newton observations appears to be 
a rather common feature in this type of AGN, is associ- 
ated with optically thick gas. This assumption needs to 
be checked in the future, at present, to our knowledge, it 
is observationally confirmed in only one of the objects in 
our sample, NGC 5506. 

We have then tested the existence of a correlation in 
R and in W a with respect to L. In both cases we find 
a marginal evidence of a very gradual increase, on aver- 
age, with increasing spectral slope. The trend in W a is 
qualitatively similar to that found in the larger sample 
of ASCA observations (gaussian line fit) by Lubiiiski & 
Zdziarski (2001, their Fig. 1), while the one that we find 
in R is much shallower than that found by Zdziarski et 
al. (1999) in Ginga data. Partly this discrepancy could be 
due to the existence of a correlation, of which we have 



found evidence in our sample, between the e-folding en- 
ergy Ef (spanning from about 80 keV to more than 300 
keV) and T, which we took automatically into account in 
our fits, while Zdziarski et al. (1999) adopted a constant 
value equal to 400 keV. Thus our assessment, based on 
an admittedly much smaller sample of objects, is that the 
same phenomenological description of the spectral com- 
ponents, which are relevant to this issue, does not confirm 
the existence of a strong correlation between the relative 
strength of the Compton reflection, and the slope of the 
power law, while it provides some evidence that a very 
gradual correlation, and only on the average, for this quan- 
tity as well as for the iron line equivalent width, cannot 
be altogether excluded, but need to be confirmed with a 
sample of objects much larger than could be studied with 
BeppoSAX. ' 

We conclude by commenting on three further limita- 
tions of this work. The first has to do with the narrow 
component of the line. On this matter further observa- 
tional steps, which could not be pursued with BeppoSAX 
alone, are required, to disentangle the broad (accretion 
disk) component from the "narrow" component of the line, 
and simultaneously the fraction of the reflection associ- 
ated with either of them. As recently shown by Malzac & 
Petrucci (2002) the nearly constant contribution by dis- 
tant and thick matter, combined with a power law spec- 
trum pivoting, as it varies in intensity, around energies in 
the range 1 to 10 keV, can introduce correlations between 
the reprocessing parameters and T. Their attempt to ex- 
plain in this scenario the paradoxical result obtained on 
NGC 5548 by Chiang et al. (2000), mentioned in Section 
3, is particularly suggestive. 

The second limitation concerns the ionization struc- 
ture and radiative transfer properties of the reprocessing 
material. Possibly with the exception of a fraction of line 
emission in NGC 5506 and Mrk 509, in our sample the 
adoption of a single zone, low ionization structure of the 
reflector seems fully justified. However, Done & Nayakshin 
(2001) have developed a disk model with a sharp transi- 
tion from outer regions of low ionization to inner regions 
with a thick skin, highly or fully ionized, where the re- 
flected continuum would be hard to distinguish from the 
direct continuum. In this model obviously the reflected 
fraction pertaining to the outer regions would present a 
value of R substantially less than unity. Moreover, their 
simulations of this model with the Ginga response, when 
fitted with the same model used by us, show that both 
R and W a turn out to correlate with L in a rather gen- 
tle manner, qualitatively similar to our results. It would 
be interesting to see this exercise repeated with the wider 
band BeppoSAX response. 

The third limitation has to do with the shape of the 
spectrum of the primary component. The conventional ex- 
ponential factor adopted is useful to reproduce the steep- 
ening of the PL at high energies, as required by the data, 
but clearly the value of R which obtains is very sensitive 
to this shape. When using instead models of the contin- 
uum based on anisotropic comptonization in a disk corona 
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(Haardt & Maraschi 1991, Haardt 1993) in fitting a sam- 
ple of BeppoSAX observations, Petrucci et al. (2001) have 
shown that the value of R turns out different, typically 
larger than that provided by a description like the one 
used by us. 

Thus on the whole there is still a very ample margin 
of progress to be made, before the confrontation between 
theory and observations could be accomplished in a more 
satisfactory way. 
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